Multi-resonance of energy transport and absence of heat pump in a force-driven lattice 
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Energy transport control in low dimensional nano-scale systems has attracted much attention in 
recent years. In this paper, we investigate the energy transport properties of Frenkel-Kontorova 
lattice subject to a periodic driving force, in particular, the resonance behavior of the energy cur- 
rent by varying the external driving frequency. It is discovered that in certain parameter ranges, 
multiple resonance peaks, instead of a single resonance, emerge. By comparing the nonlinear lattice 
model with a harmonic chain, we unravel the underlying physical mechanism for such resonance phe- 
nomenon. Other parameter dependencies of the resonance behavior are examined as well. Finally, 
we demonstrate that heat pumping is actually absent in this force-driven model. 
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I. INTRODUCTION 



In the last decade much attention has been paid to 
heat conduction on the nanoscale. Numerous studies in 
the quest of manipulating heat have brought substan- 
tial progresses in the area of Phononics, the science and 
engineering of phonons [l|. Theoretical models of var- 
ious thermal devices such as thermal rectifier 0], ther- 
mal transistor [H, thermal logic gates [ij and thermal 
memory [5[ have been proposed to control phonon-based 
thermal transport. In addition, some experimental works 
have been carried out. For instance, solid-state thermal 
diodes based on asymmetric nanotube [6|, asymmetric 
cobalt oxides Q , the nanotube phonon waveguide @ and 
solid thermal memoryQ have been realized experimen- 
tally. In the works mentioned above, heat always flows 
from high temperature region to low temperature one. 
This is in accordance with the second law of thermody- 
namics. 

However, just like its macroscopic counterpart, a heat 
pump which directs heat against thermal bias by using 
external force exists on the molecular level. Recent stud- 
ies have suggested several such kind of models based on 
different mechanisms. Li et al. proposed a heat ratchet 
to direct heat flux from one bath to another in a non- 
linear lattice, which periodically adjusts two baths' tem- 
peratures while the average remains equal 0. Later, 
Ren and Li demonstrated that heat energy can be rec- 
tified between two baths of equal temperatures at any 
instant and the correlation between the baths can even 
direct energy current against thermal bias without exter- 
nal modulations [ll|. Beyond those classical models, a 
quantum heat pump consisting of a molecule connected 
to two phonon baths was proposed recently [l| , and the 
Berry phase effect induced heat pump was unraveled as 



well [13|. 

Nevertheless, some contradiction exists in the force- 
driven heat pump. Marathe et al showed that under peri- 
odic force driving, two coupled oscillators connected with 
thermal baths fail to function as a heat pump [l4[ • Later 
Ai et al claimed the heat pumping appeared in Frenkel- 
Kontorova (FK) chain under the influence of a periodic 
driving force [15j. In addition, they observed a ther- 
mal resonance phenomenon that the heat flux attained a 
maximum value at a particular driving frequency, which 
is similar to the previously found resonance induced by 
driving bath temperatures (TlT |. However, other than the 
temperature driven clear physical mechanism of 

the driving-force-induced resonance is still unavailable. 

In this paper, we investigate the force-driven FK chain 
as a typical nonlinear lattice model and analyze ther- 
mal properties which are frequency dependent. We dis- 
cover that multiple resonance peaks, instead of a sin- 
gle peak, emerge in certain parameter ranges. The ori- 
gin of this phenomenon in fact relies on the eigenfre- 
quencies of this system. The paper is organized as fol- 
lows. First of all, we will introduce the FK model and 
briefly show the crossover from a single resonance to mul- 
tiple resonances. To look into the underlying mecha- 
nism of the multi-resonance phenomenon, we invoke a 
harmonic model through which analytic expressions of 
various quantities are obtained. Then the resonance be- 
havior of the FK model is discussed in detail. We will 
show that as far as the resonance property is concerned, 
there is much similarity between the FK model and the 
harmonic model. The effect of different system param- 
eters on the resonance phenomenon is explored as well. 
In the end, we shall briefly demonstrate that the force- 
driven model fails to perform as a heat pump. 



II. MODEL AND RESULTS 
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We start with a FK chain coupled to Langevin heat 
baths at two ends. A periodic force is applied to the left 
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f(t)=A sin(Ot) 



FIG. 1: (Color online) Schematic of a one-dimensional FK 
lattice coupled to two Langevin baths. The left end is at- 
tached to a periodic driving force f(t). 



end (See Fig. [I}. The Hamiltonian of the system is given 
as: 
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where f(t) — Aq sin(f2t) is the driving force with am- 
plitude Aq and frequency ft imposed on particle a. Xi 
denotes the displacement of the ith particle and pi the 
corresponding momentum, k is the spring constant, V 
the strength of the on-site potential and a = 1 the spac- 
ing between adjacent particles. We set a = 1 in follows 
unless otherwise stated. Under fixed boundary condition, 
the equation of motion (EOM) of the coupled particles 
(i = 2, • • • N - 1) reads: 
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The two noise terms are white noise with (T)i(t)rjj(t')) = 
2~fkBTiSijS(t — t'), = L,R), where 7 is the friction 
coefficient, ks is the Boltzmann's constant and T^ija de- 
notes the temperature of L(R) reservoir. Without loss of 
generality, we set m = 1, ks — 1 and integrate the EOM 
by the symplectic Velocity Verlet algorithm. The simu- 
lation are performed long enough (of order 10 8 ) with a 
time step of 0.005 to allow the system to reach a non- 
equilibrium steady state. The local energy current is 
defined as Ii{t) = (—kxi(xi — Xi-i)), where (• • • ) de- 
notes the ensemble average. After the transient time, 
the transport approaches an oscillatory steady state and 
the periodic average of local current is equivalent at each 
site i, reading 
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(3) 



There is no doubt that the thermal resonance exists for 
the energy current as a function of the driving frequency. 




FIG. 2: (Color online) Energy current J vs the driving fre- 
quency in the force-driven FK model. Parameters are 
T L = 0.4, Tr = 0.6, k = 1,V = 5,7 = 1, and N = 64. J 
denotes Ji (i > a = 1). The left arrow means the upper curve 
corresponds to the left scale while the right arrow means the 
lower curve corresponds to the right scale. At large Q, the 
response of the system cannot catch up with the fast driving 
as if there is no driving. At small fi, the system reduces to 
a quasi-static counterpart as if there is no driving as well. 
Therefore, at these two regions, J < since Tl < Tr. 



However, surprisingly, we observe that multiple reso- 
nance peaks emerge when the amplitude Aq is increased, 
as shown in Fig. [5J To analyze this multi-resonance phe- 
nomenon further, we invoke the harmonic chain model, 
of which clear analytic results are possible. 



Analytic results for force-driven harmonic 
lattice 



Harmonic potential is a second order approximation 
of realistic potential around its minimum. Under har- 
monic approximation, phonons will not interact with 
each other. Starting from this noninteracting picture, 
interactions within phonon gas can be introduced, by in- 
cluding higher order terms in the potential expansion. 
Compared to non-linear models such as the FK lattice, a 
harmonic model has the advantage of having an exact so- 
lution and it gives many satisfactory explanations about 
thermal properties of the crystal. 

The linear harmonic model considered in this section 
is obtained by replacing the nonlinear on-site potential 
in Eq. ([T]) by ^k a xf, where k a is the force constant of 
the on-site harmonic potential. As we shall see, under 
certain conditions the effect of the driving force can be 
separated from that of the thermal reservoirs. The EOM 
in a compact matrix form is: 



Ml = -TX + r](t) + F(t), 



(4) 



where M is the mass matrix set as identity matrix and 
$ the real symmetric force matrix, with $y = (2k + 
k D )Sij - kdt.j+i - kSi t j-i. Tij = jSij(Sa + Sin) denotes 
the coupling to reservoirs. X = (x%,X2, ■ ■ ■ ,xn) t is the 
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vector of displacements and r)(t) — (r)L,0, ■■■ ,0,r?fl) T 
depicts the thermal noise in reservoirs. F(t) is a column 
vector with elements /j(t) = Si a Ao sin(fii), denoting the 
driving force which acts on the ath particle. Defining the 
Fourier transform of quantity A and its inverse: 

-i />oo />oo 

A(u) = — e luJt A(t)dt, A(t) = / e- iujt A(uj)dw, 

2n J-oo J-oo 

and applying them to Eq. (j4]), we obtain 

X(t) = X s (t)+X d (t), (5) 

/oc 
-CO 

/OC 
duje~ lut G{u:)F(uj), 
-oo 

where G(oj) = («& — w 2 M — iujT)' 1 is the phonon Green's 
function. It is clear that the displacement vector X(t) 
is a superposition of two contributions: X s , as the ef- 
fect of the stochastic heat bath, and X d , as the con- 
sequence of the driving force. For our periodic driving 
case, it is easy to get X d {t) = -lm[G(n)fe-' int ], with 
fi = diaAo. Similarly, the velocity has the same decom- 
position X = X s + X d , with X d = lm[iflG(n)fe- int }. 
Straightforwardly, following the definition, the local en- 
ergy current can be similarly expanded as well: 

im = im+i d {t) (6) 

Jf(t) - -k(xi(xi - xU)) 

i d (t) = -kxfci-xti). 

With Ii(t) written in this form, it is obviously that the 
contributions of noise term and driving force to the en- 
ergy current are independent from each other. We should 
note this independence comes from the fact that the driv- 
ing force is not statistically correlated to the heat baths. 
Since in the present case, the driving force is determin- 
istic, we have dropped the notation of ensemble average 
(•••) in 

If(t) is just the conventional steady-state heat flux, 
which is proportional to the temperature difference and 
is given by If = ^^-^ j dujw 2 \G lN (w)\ 2 
Thus, the contribution of stochastic heat bath to the 
oscillatory steady state energy current is just Jf = 

fi/(27r) /q" Ifdt = If. In all the following discussions 
and simulations, we will set AT = Tl — Tr = un- 
less otherwise stated, since a finite temperature differ- 
ence mainly causes an additive shift of the energy current 
curve. 

To obtain the expression of Jf, we just need to calcu- 
late the product of displacement and velocity contributed 
from the driving force by applying Eq. ([3]). Denoting 
K(Q) = Re[G(fi)] and X(fi) = Im[G(fi)], we have 

x f± d = Aln[TZ ia l ja sin 2 (tit) -l ia n ja cos 2 (m) 

+ ^TZia7l ja sin(2Qi) - ^l ia l ja sin(2fii)]. (7) 



Since the driving force is periodic such that the steady 
state energy current is oscillatory, we take the periodic 
average of the physical quantities, and the final expres- 
sion of driving-force-contributed energy current is: 

kA 2 il 

Jf = -^-[Hi- la (n)Xi a (ty -zi_i a (fi)rcia(n)] 

kA 2 fl 

= -^Im[G!_ lQ (fi)G, a (fi)]. (8) 

Therefore, when T x = T 2 , J t = Jf + J d = J d . In fact, 
jf has two values {% ^ a and i > a) and is constant on 
each side of the driven site a. A direct proof is detailed 
in the Appendix. 

From the expression of Eq. ©, we should be aware 
that the energy flux possesses the denominator D(J7) = 
| det[Z(fT)]| 2 with Z(Q) = * - fl 2 M - iOT. Therefore, 
when D (f2) approaches its minimums under certain driv- 
ing frequencies, the energy flux will exhibit its maximum 
values. As a consequence, resonance emerges. Let us 
denote Pjv(f2) = det(4>Ar — £1 2 M) to be the character- 
istic polynomial of the N x N force matrix 3? at with N 
particles. It can be shown that, 

det[Z(fi)] = P N (Cl) - 7 2 fi 2 P A r_ 2 (f7) - 2i-yQP N _ 1 {Q), 

(9) 

where Pjv_i(r2) is the characteristic polynomial of the 
(N— 1) x (A — 1) force matrix <&n-i with the first row and 
column or the last row and column taken out from 3? at. 
Pn-2{Q) is the characteristic polynomial of the (A — 2) x 
(N — 2) force matrix &n-2 with both the first and last 
rows and columns taken out from 3? at. Therefore, the 
denominator D(Q) is given by 

D(p) = [p N (si) - ~f 2 n 2 p N _ 2 (n)} 2 + Aj 2 n 2 p 2 

(10) 

Apparently resonant frequencies correspond to those val- 
ues of f2 which make D(f2) locally minimized. It is dif- 
ficult to get the explicit exact solutions for those locally 
minimized solutions. Nevertheless, much simpler and in- 
spiring results can still be obtained if we consider the 
limiting case with either small friction coefficient or suf- 
ficiently large friction. For small friction, we can just set 
7 = in Eq. (TO} and obtain D(Q) = P#(0). Then, 
the resonant frequencies of the energy current are just 
the N eigenfrcquencies of the force matrix <J>at . For large 
friction, by keeping only the highest order terms of 7 in 
Eq. flTUJl, we get PJ(fi) ~ P A v_ 2 (fi). In this limiting case, 
the resonant frequencies correspond to the N — 2 eigen- 
frequencies of $at_2- When the friction is in between 
the two cases, we have to minimize D(fl) to get those 
resonant f2, which should be a gradual shift between the 
eigenfrequencies of &n-2 and $>n. 

Using Eq. l[8]). we calculate energy current vs driven 
frequencies, as shown in Fig. [3] It is not surprising 
to observe the multiple-peak resonance behavior in the 
harmonic case, since whenever the driven frequency ap- 
proaches one of the system's eigenfrequencies, resonance 
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FIG. 3: (Color online) Steady state energy current versus 
driving frequency in the harmonic lattice, k = 1, k = 0, 
A = 3.2, 7 = 1, AT = T L - T R = 0, N = 32 and a = 1. 



will occur. The number of the eigenfrequencies is mainly 
determined by the system size. 

We can also define other quantities such as the rate of 
heat released from the two heat baths and rate of work 
done by the driving force. Their definitions and analyti- 
cal expressions are given below: 

q L = ( Xl ( VL - yxx)) = (±t( VL - 7 iJ)> - 1 (x d l )\ (11) 

where stands for the rate of heat released from the 
left bath. The first term is contributed by heat baths 
and is zero when AT = 0, while the second term comes 
from the driving force. Therefore the periodic average at 
AT = gives 

O r 2 */^ -yA 2 Q 2 

Ql = ^L dt * L = ® l + ^ = - 2 ^r-\ G ^)\ 2 . 

(12) 

Similarly, the expression for Qr, the average rate of heat 
released from the right bath, is given as follows 



Qr — Qr + Qr 



t4P 

2 
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(13) 



The rate of work done by the driving force is defined as 
w(t) = (f(t)x a ) = f(t)i d a (14) 
where f(t) = Aos'm(ilt). Its periodical average gives 

h) 



O r 2 ^l a A 2 Q 

w = — dtw(t) = -^-i aa (n). (15) 

2tt J a I 



Just like the energy flux, W, Ql and Qr also show multi- 
resonance behavior jl7| . By the continuity equation, the 
energy current flowing out and into a particle must can- 
cel each other when the system reaches its steady state, 
since the local energy density does not vary with time 
at steady state. Therefore, we immediately obtain the 
relation W + Ql = Ji = —Qr (i> at), which obeys the 
energy conservation. For the harmonic model, a direct 
proof is given in the Appendix. We should also be aware 
that this relation also holds in the FK model and this is 
verified by numerical results. 



FIG. 4: (Color online) Comparison of the FK case and har- 
monic case with force amplitude Ao = 3.2 in a randomly se- 
lected frequency region. Other parameters are k = 1, k a = 0, 
V = 5, 7 = 1, AT = 0, N = 64 and a = 1. 



B. Multiple resonances of energy transport in 
force-driven FK model 



The main qualitative difference between the multiple- 
resonance curves in Fig. [2] and [3] lies in the fact that the 
height of the peak for the FK case is smaller than that 
for the harmonic case, due to the nonlinear on-site po- 
tential. This is more obvious for low and high frequency 
regime. More importantly, we notice that the positions 
of multiple peaks of the two models seem to be close 
to each other, as illustrated in Fig. This suggests 
that the FK model and harmonic model share a similar 
resonance mechanism. The harmonic model discussed 
in Sec. Ill Al is able to shed lights on the multi-resonance 
mechanism in a force-driven nonlinear lattice. It is worth- 
while to point out that when sampling frequency is not 
dense enough, one may get a single peak, which actually 
corresponds to the lower envelope of the multi-resonance 
curve. When Ao decreases, the system seems to have a 
single peak even though we increase the density of the 
sampling points, as shown in Fig. [5] The mechanism for 
the single peak should still be explained by the fact that 
the driving frequency is resonant with the system's char- 
acteristic frequencies. However, the non-linear potential 
smooths out the peaks of multi-resonances. In this way, 
we can only observe the lower envelope so that only one 
peak is observable. Thus the single-resonance curve does 
not completely reflect the intrinsic transport property of 
the FK model. 

Let us examine how the parameters of the system will 
affect the multi-resonance curve. We have already seen 
the case of the appearance of multiple peaks when Aq in- 
creases, as shown in Fig. [2] When Ao becomes large, the 
kinetic energy of the particle gained from driving force 
will be much larger than the height of the on-site po- 
tential. In this way, the mask effect of nonlinear on- 
site potential becomes negligible and the FK model ef- 
fectively reduces to a harmonic model without on-site 
potential. This is also demonstrated in Fig. [S] with 
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FIG. 5: (Color online) Heat current vs driven frequency for 
large force amplitude Ao = 10. (a) Low frequency regime, 
(b) Moderate frequency regime, (c) High frequency regime. 
Other parameters are the same as in Fig. 0] 
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Aq = 10. We arbitrarily choose three typical regimes 
of driving frequency: small, moderate and large fi. In 
all three regimes, the FK model is very close to the har- 
monic one. Thus, as long as A becomes large enough, 
the mask effect of nonlinearity will fade away and the 
intrinsic multiple resonant peaks will emerge. 

We are also interested in the effect of temperature. 
Considering that T (defined as T = (T L +T R )/2) will 
take the role of thermal excitation which mainly affects 
the kinetic energy of particles, we expect a large value 
of T will have a similar effect with A$ on the resonance 
behaviors. As shown in Fig. [5] with T — 0.5 and Aq = 1, 
there is no multiple resonances behavior. While for the 
increased temperature case T = 2.0, there is clearly a 
peak in a small f2 range, as shown in Fig. H^a), which 
indicates the multiple peaks resonance behavior in the 
whole frequency range. At low temperature with large 
amplitude, this is also the illustrated in Fig.[5][b). 

When T and Aq are both small, compared to the nonlin- 
ear on-site potential V, particles are confined near their 
equilibrium positions, within the valley of non-linear on- 
site potential. By Taylor expansion of the potential to 
the second order, we see that the non-linear potential 
reduces to a harmonic one, so that the harmonic ap- 
proximation works. In this case, the lower bound of the 
phonon band will be raised by yV. The phonon band 
in this case is therefore moved from < Q < \/4k to 
W < < \JV + 4fc [l8j]. (Accordingly, the phonon 
band of the harmonic model is shifted by the on-site po- 
tential k a to \fk < VL < \/k Q + 4k). Energy transport 
is forbidden outside the frequency region. It is thus ex- 
pected to see a shift of the energy flux curve to the right 
while the multi-resonance are still observed, as shown in 
Fig. He). 

The above discussion reveals that there are mainly 
three dynamic regimes: 

(1) When k B T + A a/(2ir) > V/(2ir) 2 (a is the lat- 
tice constant as defined previously), the FK model will 
approach a harmonic model without an on-site potential. 



FIG. 6: (Color online) Temperature effect on comparison of 
the FK case and harmonic case, (a) T = 2.0, Aq = 1. (b) T = 
0.01, A = 6.4. (c) T = 0.01, Ao = 0.1. Other parameters 
are the same as in Fig. 2] except in (c): k = 10, k — 15 and 
V = 15. 



Multi-resonances are observable. 

(2) For the opposite case, in which fc^T + A a/(2Tr) <C 
V/(2ir) 2 , the FK model reduces to a harmonic model 
with pure harmonic on-site potential of strength k a = V . 
Multi-resonances are still observable. 

(3) When T + A a/{2n) is comparable with V/(2tt) 2 , 
the multiple peaks will be smoothed out by the effect of 
nonlinear on-site potentials and only single resonant peak 
is observable. In the intermediate regimes, there is the 
crossover from single peak to multiple peaks, of which the 
resonant magnitudes are smaller compared with those in 
the harmonic model. 

In addition, the heat released from the baths and the 
work done by the force in FK model also show parameter- 
dependent multiple resonances [17] , which can be ex- 
plained by the same mechanism. 



III. ABSENCE OF HEAT PUMPING IN A 
FORCE-DRIVEN LATTICE 

In all simulations so far, we only consider the energy 
flux flowing into the right reservoir as a whole. In fact 
it is a sum of two parts: the rate of heat released from 
the left reservoir into the system and the rate of work 
done by f(t), i.e, J — —Qr = Ql + W. Now let's con- 
sider a typical result of energy transport in a force-driven 
FK lattice with Tl < Tr, where these two contributions 
are separated, as shown in Fig. [71 Notably, though J 
is positive (from L to R) in the resonance region, Ql is 
negative in the full range. This indicates heat will always 
flow into the left reservoir from the chain, whatever the 
driving frequency is. The corresponding energy flow di- 
agram is depicted in Fig. [S^a). This actually shows that 
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FIG. 7: (Color online) —Ql, W, J vs the driving frequency 
in force-driven FK model. Parameters are the same as those 
for lower single peak curve with Ao = 1, Tl < Tr in Fig. [2] 
Ql is negative at all frequencies, showing a negative result 
for pumping effect. 




is always negative since Tl < Tr. And the second driv- 
ing contribution, which docs not depend on the temper- 
ature, always contributes a negative value as well. Thus, 
Ql will still be negative, which indicates the low tem- 
perature bath always absorbs rather than releases heat, 
even we may change the position of the driving force. 
The situation just corresponds the energy flow diagram 
in Fig.lUJa). 

Clearly as long as the effects of the driving force and 
the baths are independent, Eq. (jlip holds for harmonic 
chain. Thus for a harmonic system, in order to get a 
heat pump, a necessary condition is that the driving force 
is statistically correlated to the bath, such that these 
two contributions could be synergetic. While in the FK 
model, effects from the deterministic driving force and 
stochastic baths are not separable as in Eq. (fTTj) because 
of the nonlinear on-site potential. However, the non- 
linearity can not yet make the two contributions syner- 
getic. And still, our numerical simulations show a nega- 
tive result for a force-driven FK lattice acting as a heat 
pump. We thus speculate that for both the force-driven 
harmonic and nonlinear lattice, the correlation may be 
the key ingredient for the presence of heat pump ac- 
tion, similar to the entanglement in quantum thermal 
baths [l9l | , the off-equilibrium nonthermal reservoirs (20j 
and the bath noise correlation in temperature-driven case 
11], for "low temperature to high temperature" thermal 
transports. 




FIG. 8: (Color online) (a) A typical energy flow diagram in 
force-driven FK model with Tl < Tr. (b) The situation when 
a model performs as a heat pump with Tl < Tr. The arrow 
denotes the direction of the energy flow. 



the FK system fails to function as a heat pump since nor- 
mally, a pump should have an energy flow diagram as in 
Fig. [EJb), where the heat is absorbed from the low tem- 
perature bath and released to the high temperature one. 
We have tested a broad range of parameters and even 
in the regimes with multiple resonances. In all cases, the 
energy never flows as in Fig. [8jb) . Therefore we conclude 
that there is no heat pumping action in a force-driven FK 
lattice. 



This finding is consistent with Ref. wherein it is 
found that two coupled oscillators of either harmonic or 
FPU-like interaction would not act as a heat pump un- 
der external driving forces. Ref. [IU gets the opposite 
conclusion, because they overlook the direction oi Ql 
and claim the model can act as a heat pump whenever J 
flows to the right bath of a higher temperature. 

To understand the underling mechanism for the ab- 
sence of heat pumping in a force-driven lattice system, we 
may consider Eq. (|lip . The first stochastic contribution 



IV. CONCLUSION 

To summarize, we have studied the energy trans- 
port control in a force-driven one dimensional nonlin- 
ear lattice-the Frenkel-Kontorova model. We have found 
multiple thermal resonances as a function of the driven 
frequency. Although not exactly the same, the resonant 
frequencies are closely related to the eigenfrequencies of 
the force matrix, and the number of resonant peaks is 
bounded above by the system size N. Moreover, since 
the onsite nonlinear potential tends to decrease the mag- 
nitude of energy current and smoothes out the multiple 
peaks, the multi-resonance phenomenon is only observ- 
able in certain parameter ranges and the crossover from 
multiple resonances to single resonance will occur. Fi- 
nally, by following a rigorous definition of heat pump, we 
clarify a previous contradiction and conclude that heat 
pumping effect is absent in force-driven lattices. 

In this paper, we focus on a one-dimensional chain. 
However the analysis and results can be generalized to 
high dimensional system with arbitrary topology, like 
polymer networks, proteins [2lfl. and the three dimen- 
sional random elastic network [22|], of which the eigen- 
spectra are much more complicated and non-trivial. The 
investigation of the resonances in such systems will give 
us more flexible methods for mechanical control of energy 
transport in real applications. 
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APPENDIX 

Intuitively, from the energy conservation point of view, 
the following relations should hold 



-Qr — Ja+l — ■ ■ ■ — Jn, 
Ql = J\ = ■ ■ ■ = Ja, 

W = -Q l -Qr. 



(16) 
(17) 
(18) 



However, it is highly nontrivial to give direct proofs for 
arbitrary nonlinear force-driven lattices. In the follow- 
ing, we will provide a direct proof of the above equations 
in force-driven harmonic chains by using the analytical 
forms of Ql,r, W and J,. 

Let us first define the dimensionless matrix Z(fi) = 
Z(f2)/fc, such that = 5i.j {c — ib5n — ibS^) — 5i,j+\ — 
6ij-i, where b — jCl/k and c = 2 — mtt 2 /k are two 
dimensionless parameters. Here we set k a = without 
loss of generality. Considering G = Z _1 = Z _1 /fc, we 
need a formula for the inversion of the matrix Z, which 
is given in its recurrence form [23l |: 



?i-l<Pj+l 

9j-lfa+l 
~0n 



if i < j, 
if i > j, 



(19) 
(20) 



where 9i and fa are two sequences, defined by the recur- 
rence relations: 



— Zadi-i — 0i-2, 



(21) 
(22) 



with the initial condition 8q = 1, 0\ = Zn, 4>n+i = 1 
and 4>n — Znn- Note 9n is actually the determinant of 
the matrix Z, which will be denoted as d in follows. Now, 
to prove Eq. (fT6|) , let us recall the expression for Qr and 
Ji (see Eq. © and |J3J|) and use Eq. O and O- We 
then get (set Aq = 1 for clearness): 



7 r> 2 



*- ! T |G, ' a|2 -«¥ l 

fcfi 



(23) 



2k\d\ 2 
illm[(j)*fa + i 
2k\d\ 2 



lm[(0 a -ifa)*(9 a -ifa+i)] 



(24) 



where i > a. Hence, to prove the above two equa- 
tions are equivalent with each other, we need to show 
lm[0*0i + i] = yfl/k for all particles to the right of the 
driving force. We do induction on the index /. For I = N, 
the relation is clearly true, which can be taken as the 
base case. Then suppose the relation is true for all I > i, 
we need to show it is also true for I = i — 1 case. Use 
the recurrence equation Eq. (|2ip . we have Im^*^^] = 
Im[(cfa*-$ +1 )fa] = -Im[$ +1 fa] = Im^i+i] = -yQ/k. 
Thus, we have finished the proof that Eq. {jTHJ) holds for 
all particles to the right of the force. A similar proof can 
be given to show Eq. (TTT)) holds for all particles to the 
left of the driving force. 

Before proceeding to prove Eq. (fT8|) . we need more 
notations. Denote the n x n matrix M^") with elements 
M^ n \j = cSi t j — ^ij+i — 5i,j-i and its determinant K n = 
det[M^], with boundary K-i = and K = 1. The 
explicit expression of the determinant is given by K n —\ = 
(Vi - V2 )l(Vi - 2/2), where yi, y% is the solution of the 
quadratic equation y 2 — cy + 1 = 0. We observe that 0, 
is in fact the determinant of the sub-matrix in Z starting 
from the 1st row and column to the ith row and column. 
Similarly fa is the determinant of the sub-matrix in Z 
starting from the ith row and column to the Nth row and 
column. Now we may expand 9i , fa and the determinant 
of Z in terms of K n 



= K N - b 2 K N - 2 - 2ibK 
= Ki - ibKi-u 
= Kn-i+i — ibKN--, 



\ -l, (25) 
(26) 
(27) 

Now we are ready to prove Eq. fTS")) . Let us assume 
Eq. (TTg|) holds and substitute in the expressions of W, 
Q L and Q R (see Eq. ([T5 ]) . (p ]) and {IS ]) ). Then Eq. iflS]) 
becomes l aa = jn{\G la \ 2 + \G Na \ 2 )- With b = jtt/k 
and G = Z^ 1 /k, also using Eq. (fT9])~ ([22]). we get 

lm(e a ^<f> a+1 d*) = b{\e a ^\ 2 + |0 Q+ i| 2 ). (28) 

Express the left hand side (LHS) and right hand side 
(RHS) of Eq. (H in terms of K n (see Eq.(25)~(27)), 
and after some algebraic work, it can be shown 

LHS =b 3 (K a - 2 K N - 2 KN-a + Ku-\Kn-2K N-a-l 

— ^Ka-2KN-\K N-a-l) 

+ b{2K a ^iK N -\K N - a — K a -2KNK N - a 

- K^KnKn.^), (29) 
RHS =b\K 2 a _ 2 + K N _ a _ x ) + b{K 2 a _ x + K 2 N _ a ), (30) 



Now in order to show Eq. ()T8|) holds, it is sufficient to 
prove LHS=RHS to complete the proof. Denote the coef- 
ficient of b 3 (b) in LHS by L 3 {L x ) and that of RHS by R 3 
(Ri). By using the formula K n -i = (y{ 1 - y2)/(yi - 2/2) 
and the fact yiy2 = 1, we get 



2a-l 2a-l . 2N-2a + l 2N-2a+l 

3 nil -V2 +Vl V2 



(yi - JJ2 

I 2a — A 'In- s, , 1 

-\Vi - 2/2 ) - [Vi 
- 4(yi - 2/2)] = Rs- 



2N-2a-l 



2N-2a-l 
2/2 



(31) 
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Note that by a relabeling of a — > a + 1 and N — > ./V + 2 K n _2, we can show L3 — > Li and i?3 — > Therefore, 
in L 3 and i? 3 , and utilizing the relation K n = cK n _\ — LHS=RHS and the proof is completed. 
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